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Exchange interactions in a- and 7-Fe are investigated within an ab-initio spin spiral approach. 
We have performed total energy calculations for different magnetic structures as a function of lattice 
distortions, related with various cell volumes and the Bain tetragonal deformations. The effective 
exchange parameters in 7-Fe are very sensitive to the lattice distortions, leading to the ferromagnetic 
ground state for the tetragonal deformation or increase of the volume cell. At the same time, the 
magnetic-structure-independent part of the total energy changes very slowly with the tetragonal 
deformations. The computational results demonstrate a strong mutual dependence of crystal and 
magnetic structures in Fe and explain the observable "anti-Invar" behavior of thermal expansion 
coefficient in 7-Fe. 
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I. INTRODUCTION 

Iron-based alloys are still among the most important 
industrial materials. The thermodynamic properties and 
mechanism of phase transformations in these materials 
have been discussed intensively last yearsi. Nevertheless, 
the fundamental properties of iron have not been com- 
pletely understood up to now. Main difficulties are re- 
lated with a non-trivial combination of the itinerant and 
localized behavior and correlation effects of 3d-electrons 
determining electronic, magnetic and structural proper- 
ties of iron^r— . It is commonly accepted now that mag- 
netic degrees of freedom play a crucial role in the stabil- 
ity of different iron phases^— which makes the situation 
even more complicated. Interplay between magnetic and 
lattice degrees of freedom in different crystallographic 
phases of iron remains still unresolved problem. 

One of the most complicated example related 7-phase 
of iron with highly frustrated magnetic structure. There 
are many magnetic configurations of 7-Fe with almost 
the same total energies and the ground state is crucially 
depends on the value of lattice parameters^— . The sen- 
sitivity to dilatation has been studied in detail by many 
groups^— in the context of a so called moment- volume 
instability^ 7 -. At the same time, the energy dependence 
on the tetragonal deformation which is closely related 
with the Bain deformation path of a — 7 phase transfor- 
mation also deserves a serious attention. We discussed 
this issue in our previous worki^ and found that the 
transition of 7-Fe to the ferromagnetic state can trigger 
the martensitic transformation without noticeable energy 



barriers. In more detail, the effect of tetragonal deforma- 
tions on magnetism and vice versa was discussed in rela- 
tion with the Invar behavior observed in Fi-Ni alloy o 17 ' 19 . 
A magnetoelastic spin-lattice coupling plays also an im- 
portant role in structural phase transitions in 
and Cr-based alloys^, as well as in the magnetic shape- 
memory alloy Ni2MnGa22. A soft-mode phonon behav- 
ior, as a precursor of the 7 — » a transformation, was 
recently observed in Fe-Ni alloys^. 

In contrast with Fe-Ni alloys, the equilibrium 7-phase 
in pure Fe exists only at high temperatures T > 1200 K 
where thermal fluctuations are very strong and magnetic 
moments are disordered. Observation of the so-called 
"anti-Invar" behavior of 7-Fe^ can be related with the 
fact that the spin-lattice coupling is strong enough to 
affects the thermodynamic properties up to very high 
temperatures. 

In this paper we investigate quantitatively a variation 
of the exchange parameters in a- and 7-Fe as functions 
of tetragonal Bain-deformations and dilatation. Whereas 
the sensitivity of the exchange parameters to the di- 
latation has been studied previously^— an information 
about the tetragonal deformations have been missing un- 
til now. Based on the calculated magnetic exchange data 
we discuss the origin of the anti-Invar behavior of 7-Fe. 



II. COMPUTATIONAL APPROACH 

The standard approach to study magnetic properties of 
itinerant-electron transition-metal systems related with 
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the mapping of density functional total energies on the 
effective classical Heisenberg model: 

i<j 

were e,; is the unit vector in direction of the magnetic 
moment at site a' 25 ' 26 . In this notation the value of on- 
site atomic magnetic moments Mi is included into the 
exchange parameters Ji j. Therefore the total energy of 
the system is a sum of a magnetic-structure independent 
contribution Eq and the "Heiscnbcrg-cxchange" part: 
E = Eq + H ex , where Eq is a function of deformations 
and the magnitude of local moments: Ep(fl, c/a, M). A 
similar decomposition was used earlier in Ref. 1271 . One 
should stress that Eq is dependent on the values of mag- 
netic moments Mi and is therefore essentially different 
from the energy of a non-spin-polarized state, which at- 
tribute is zero all magnetic moments: Mi = 0. 

There are two main approaches to the mapping onto 
magnetic Hamiltonian. An analytical scheme is based on 
the use of so-called "magnetic force theorem" 25 ' 26 , when 
the exchange interactions are obtained from variations 
of the total energy with respect to infinitesimal devia- 
tions of the magnetic moments from a collinear state. In 
this paper we use more accurate numerical method based 
on the density functional calculations of the spin spiral 
magnetic structures where the neighboring magnetic mo- 
ments are rotated relative to each other by a finite angle 
(for review, see Ref. |28| ). This scheme includes a spin- 
and charge-density relaxation for large moment fluctua- 
tions. The energy per atom of the spin spiral with the 
wave vector Q can be presented as: 

= £ -^Z„J n exp(iQ-R„), (2) 

n 

where N is the number of magnetic atoms, Z n is the num- 
ber of the n-th nearest neighbor atoms, Eq is a magnetic- 
structure-indcpcndcnt contribution to the total energy of 
the system, Rjj is the vector connecting sites i and j, 
n labels the coordination shell. The exchange parame- 
ters J n can be found from Eq. © by using the discrete 
Fourier transformation: 

J„ = -i^ J E(Q fe )exp(iQ fc -R n ), (3) 

k 

where the summation runs over a regular Q-vcctor mesh 
in Brillouin zone with the total number of points K. As 
follows from Eq.Q the value Eq is the average value of 
spin spiral energies over all Qfc, 

Eo = ^Y, E (^y ( 4 ) 

k 

In principle, one can find the dependence of total energy 
on magnitude of M within a constrained moment spin- 
spiral calculations, but this lay beyond a scope of present 



paper. A parameter of total exchange energy 

Jo = /J Z n J n (5) 

n 

characterizes a ferromagnetic contribution to the total 
energy. Note that the decomposition of the total energy 
used in Ref. [27] differs from that used in this work by a 
shift by Jo- 

In general, the exchange parameters found from the 
planar spin spiral calculations and from the magnetic 
force theorem arc different and only the value of a spin 
stiffness constant should be the same 2 ^. Note that 
parameters obtained by the use of infinitesimal spin 
deviation o 25 ' 26 give a correct description of a magnon 
spectra, while parameters found from a direct calcula- 
tion of the spin spiral total energies are supposed to 
be more accurate for descriptions of thermodynamic 
properties 2 ^. The difference of J n obtained within these 
two approaches characterizes a non-Hciscnberg charac- 
ter of magnetic interactions which is expected for itin- 
erant magnets such as iror* 3 ^. Another manifestation of 
the non-Heisenberg behavior related with the fact that 
the magnitude of the magnetic moments dependent on 
the spin spiral wave vector Q. Therefore, the values J n 
obtained in the framework of spin spiral approach are 
considered as effective exchange parameters. 

The total energy calculations of Fe with spin spirals 
magnetic structure is performed using VASP (Vienna Ab- 
initio Simulation Package )2ir— with first-principle pseu- 
dopotcntials constructed by the projected augmented 
wave method (PAW)21 Following an experience on non- 
collincar magnetic investigation^ we employed the gen- 
eralized gradient approximation (GGA) for the density 
functional in a form by Perdew and Wang (1 991)25. with 
the spin-interpolation^ 6 -. The PAW potential without 
core states and with energy mesh cutoff 530 eV, and the 
uniform fc-point 12x12x12 mesh in the Monkhorst-Park 
scheme^ with 1728 fc-points are used. The calculations 
are done for a single-atom unit cell subjected by two 
types homogeneous deformations, namely, dilatation (a 
change of the volume for a fixed c/a ratio) and tetragonal 
ones (a change of c/a ratio at a fixed volume). For given 
lattice parameters, the energy set E(Qk) is calculated on 
a uniform 16x16x16 mesh and the Fourier transforma- 
tion (Eq.[3|) is used to determine the exchange parameters 

Jn- 



III. COMPUTATIONAL RESULTS 

The local magnetic moments M(Q) and total ener- 
gies for the spin spiral states E(Q), calculated for dif- 
ferent values of volume and tetragonal deformations are 
presented in Figs. Q] and [2 respectively. We show the 
results only for the symmetric directions of the wave vec- 
tor r — Z — W(U) in Brillouin zone parallel to (001) and 
(012) in lattice with cubic (tetragonal) symmetry^. 
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FIG. 1: (Color online) The magnetic moment of iron as func- 
tion of the spin spiral vector Q and the lattice deformations. 



The magnetic moments depend strongly on the spin 
spiral wave vector Q, as one can see from Fig. [TJ This 
fact confirms the non-Heisenberg character of magnetic 
interactions in iron. The magnitude of magnetic mo- 
ments gradually decrease by about 30% along T — Z di- 
rection for all considered structures except for fee iron 
at small volume fl = 11.0 A 3 . Large difference of mag- 
netic moments for fee Fc at T-point at a small volume 
{n = 11.0 A 3 ) and bigger ones (0 > 11.44 A 3 ) results 
from a well known magnetovolume instability which was 
discussed in the context of the Invar problem^. 

According to our results (Fig. [5]) the ground state of 
fee iron is spin spiral with Q varying nearby 0.5(001) (in 
2tt/c units) with volume and c/a ratio for a broad inter- 
val 10.5 < < 12.0 A 3 . The magnetic ground state of 
fee iron is a controversial issue up to now. The antifcrro- 
magnetic double layer structure (AFMD), equivalent the 
spin spiral with 0.5(001) has been discussed in a series of 
paper a n ' 48 i 49 . The later publication o 13 ' 14 i 51 show, rather 
incommensurate ground states with Q-vector depending 
on lattice parameters. Our result are in agreement with 
the recent calculations 1 ^ 4 * 3 ^! 4 ^. 

An increase of iron volume further fl > 12.0 A re- 
sults in the transition from spin spiral to ferromagnetic 
(FM) structure (Fig. [2^). The energy difference AE M be- 
tween FM and antiferromagnctic (AFM) states (or spin 
spiral structure with Q = (001)) gives a scale of the ex- 
change interaction energy which decreases monotonously 
with increasing of the volume and finally changes the sign 
near n cxp = 11.44 A 3 . This volume corresponds to an 
experimental value for precipitates of 7-Fe in Cu at low 
temperatures^. 

Our results demonstrate that the magnetic structure of 
fee iron is strongly dependent on the lattice deformations 
(Fig. [5]). This conclusion agrees well with the previous 
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FIG. 2: (Color online) Dependence of the total energy of 
fee Fe on the the spin spiral wave vector -B(Q) for differ- 
ent volumes (upper panel, c/a—1) and for different c/a ratios 
(lower panel, 0=11.44 A 3 ). The deformations c/a = \/\f2 
and c/a = 1 correspond to the bcc (a-Fe) and fee (7-Fe) 
structures, respectively. Symbols mark the results of the den- 
sity functional calculations whereas the lines correspond Eq. 
(J2j) with obtained exchange parameters J„. 



investigations of iro n 13 i 18 i 51 . In particular, the spin spiral 
ground state is changed to the ferromagnetic one within 
the tetragonal deformation region along the Bain path 
from the fee (c/a = 1) to bcc iron (c/a = l/\/2). A 
magnetic transition to the FM state and its role in the 
martensitic transformation have been discussed earlier in 
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Ref. [TH. In the opposite case, when c/a > 1, the tetrag- 
onal deformation leads to a weaker dependence of E(Q). 
The spin spiral structure represents a ground state at 
c/a f=s 1.1 and a transition to the ferromagnetic ground 
state appeared at c/a > 1.2. This magnetic transforma- 
tions are in agreement with a previously obtained phase 
diagram^. 

The results for exchange parameters J n {c/a, f2) as 
function of lattice distortions are presented in Fig. [3] 
Positive values indicate that the ferromagnetic type of 
ordering is preferable. The dependence -E(Q) determined 
by Eq. ([2]) with obtained exchange parameters J n give a 
perfect interpolation to the calculated spin spiral energies 
(lines and symbols in Fig. [2]). A striking feature of this 
curves is that the total exchange energy Jo behaves simi- 
lar to Z\ Ji [Zi corresponds to number of i-th neighbors ) 
for all deformations considered. This means that the con- 
tributions of longer-range exchange interactions (n > 1) 
are canceled out. Similar results have been obtained by 
analytical calculations of exchange parameters^ for the 
volume variation of fee iron. 

Effects of volume variation on the exchange parameters 
in fee structure is very noticeable and J\ demonstrates 
there a non-monotonous behavior (Fig. [3^). At low vol- 
umes (O < fiexp) total exchange parameter Jo become 
negative showing the tendency to antiferromagnetic-type 
coupling. For atomic volumes near f2 eX p the parameter 
Ji is close to zero and the exchange energy Jo is small 
and negative. In this case the value Jo is determined by 
all exchange parameters J n with n > 1. Therefore, com- 
putational results for Q sa f2 oxp appear to be quite sensi- 
tive to the details of the approximation used^ (e.g. the 
exchange-correlation functional, energy cut-off, number 
of k-points, etc.). Such behavior of exchange parameters 
can likely be related to a complex magnetic structure 
discussed in the experimental work by Tsunoda and co- 
workers^. 

Parameter J changes the sign at a volume which is 
just slightly above the Q cxp and grows rapidly, therefore 
the ferromagnetic order becomes more stable for higher 
volumes (Fig. [3^). The behavior of J (Fig. [2]) agree well 
with previous calculation o 15 ' 16 and reproduces the tran- 
sition from spin spiral to ferromagnetic state discussed 
above. 

In the bee Fe exchange parameters demonstrate a 
rather weak sensitivity to the volume variation and the 
nearest neighbor contribution Ji is large, positive and 
dominant in a broad interval of fl (Fig. As a re- 

sult, the ferromagnetic ground state has an essential 
preference in bee Fe in comparison with the AFM and 
non-collincar magnetic structures. Results of previous 
calculation o 26 ' 39 " — give slightly lower values Jo and Ji 
in bee Fe than obtained here but also reproduce a domi- 
nate contribution of J\ to the exchange energy. 

A dependence of exchange parameters on the tetrag- 
onal deformation c/a is presented in the Fig. U3c,d). A 
symmetry break caused by tetragonal deformations leads 
to a modification of the coordination numbers in fee or 



bee lattice. Here we neglect the rearrangement of site 
positions and assume that the set of atoms belongs to 
the same coordination shells n in fee and fct structures 
for 0.85 < c/a < 1.2 and in bee and bet structures for 
0.6 < c/a < 0.85. The curves J\{c/a) and Jo(c/a) have 
both a minima for fee and a maxima for the bee struc- 
tures. One can see that near the bec structure Jo is much 
less sensitive to the dilatation than to the tetragonal de- 
formation. Near the fee, Jo is very sensitive to both types 
of deformations. This is mainly due to sensitivity of Ji to 
deformations whereas J n for n > 1 arc almost unchanged 
with variation of lattice parameters. 

The dependence of the exchange parameters on both 
types of deformations is shown in Fig. [4] as a contour 
plot Jo(fl,c/a). One can see that the tetragonal defor- 
mation together with the increase in volume enhance sig- 
nificantly the exchange interaction energy in 7-Fc. The 
value Jo ~ 70 mcV is reached for the experimental vol- 
ume of 7-Fe Q m 12 A 3 and (c/a - 1) « 5%. 

Calculated exchange parameters J„ are presented in 
the Fig. [5] as functions of interatomic distances. The 
exchange interactions in fee iron have a very long-ranged 
behavior at the volumes 57 ^ G xp- Such a strong Fricdcl 
oscillations was already found in Ref. [l6|. This is a rea- 
son of magnetic frustrations and existence of numerous 
complex magnetic structures with low energies in the fee 
F o 12 ' 16 ' 38 ' 50 . A tetragonal deformation of the fee struc- 
ture changes dramatically the behavior of J n due to a 
sharp increase of J\ contribution which becomes a dom- 
inant one. The increase of volume acts in a similar way. 
One can see from Fig. [5] that the exchange interactions 
depends not only on interatomic distance R n but also 
very sensitive to particular values of c/a. Therefore, cor- 
rect lattice deformations should be necessarily taken into 
account explicitly for a correct description of magnetic 
structures in Fe. 



IV. DISCUSSION AND CONCLUSIONS 

We can determine the magnetic-structure independent 
contribution Eq by subtracting the Hcisenbcrg-like con- 
tribution with calculated exchange parameters from the 
total energy. In order to do this one can use the en- 
ergy of ferromagnetic state in the spin-spiral framework 
Q = from the Eq. ([2]) and the following expression: 
Eq = Efm — Jo- As was mentioned earlier, Eq essen- 
tially differs from a total energy obtained in the non-spin- 
polarizcd calculations because of implicit dependence of 
Eq on the magnetic moment M. They are equal only for 
the systems with zero magnetic moments of all atoms. 

The results for Eq are shown in Fig. [5] together with 
the total energies of FM bee and FM, AFM and AFMD 
fee states obtained by the reconstruction from Eq and 
J n . For comparison, the total energies of -Enm obtained 
from calculations by VASP are also shown. These results 
agree very well with the previous ab-initio calculations^ 9 - 
and demonstrate the dramatic difference between Eq (see 




FIG. 3: (Color online) Exchange parameters J n for n = 1, 2, 3, 4, 5 for different lattice parameters: dependence J„ on a volume 
of fee (a) and bec (b) Fe; dependence J n on (c/a) at fixed volumes Q. = 11.44 A 3 (c) and O. = 12.0 A (d), respectively. 




c/a ratio 



FIG. 4: (Color online) Dependence of the total exchange 
parameter Jo on volume fl and c/a ratio as a contour plot 
Jo (fi, c/a). 




FIG. 5: (Color online) The exchange parameter as a func- 
tion of interatomic distance to the n-th neighbour J n (R n ) for 
different c/a ratios. 



Eq.SD and E NM - 

For fee iron the magnetic-structure independent con- 
tribution Eq is rather close to the energy of AFM and 
AFMD states. For bec iron the difference between Eq 
and ground-state energy E-pm is larger but rather weakly 
volume dependent compare to fee states. At the same 
time, the energy of FM fee state shows two minima at 
low and high volumes. This behavior of fee total energy 
drastically differs from the Eq curve. The difference is 
larger for higher volumes and has entirely magnetic ori- 
gin due to increase of the exchange parameters with f2 
(Fig. [3]) . Quantitatively, the values of bulk modulus for 
fee iron obtained from the Birch-Murnaghan equation of 
stated for Eq and £afmd, -Eafm curves differ by about 
17% and 30% (161, 189, and 207 GPa, respectively). For 
the bec iron estimation of bulk modulus from Eq and 
E-pM curves give the same B w 187 GPa, in agreement 
with the experiment^. 

The situation with tetragonal deformations is quite un- 
usual. One can see in Fig. |6Jd that Eq depends on c/a 
very weakly. This means that the Heisenbcrg-likc contri- 
bution is dominant in the shear modulus C , as well as 
in the whole energy curve along the Bain path. This is 
main origin of anomalously strong coupling between the 
magnetic and lattice degrees of freedom in iron, where 
the tetragonal deformation plays a special role. The 
curve Eq(c/cl) has a minimum at c/a = 1 (fee structure) 
whereas both £fm and -Eafm have no minima at this 
point which means instability of fee phase in both mag- 
netic structures. The minima correspond to bec (FM) 
and fct (AFM, AFMD) states with c/a > 1. 
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FIG. 6: (Color online) Total energy of iron reconstructed from the Eo and J n parameters (symbols) as a function of volume 
(a) and tetragonal deformation (b) for FM and AFM magnetic structures, as well as the magnetic-structure-independent 
contribution Eo to the total energies. The lines are total energy fitting using Birch-Murnaghan equation of stated for (a) and 
splines for (b). Fitting for the FM state in fee structure (a) is done separately for each minimum. 



Our calculations reveal another unusual feature of the 
magnetic interactions in fee iron related with a growth of 
the exchange parameter J\ and, as a consequence, Jq 
with the volume increase at Q > f2 exp (see Fig. [3^,). 
This behavior corresponds to the rising branch of the 
Bethe-Slater curve J(f2) which have been used for a 
semi-quantitative interpretation of the Invar anomaly^. 
This region of volumes corresponds to observed high- 
tcmpcrature phase of 7-Fe; for further increase of inter- 
atomic distances the overlap of c?-orbitals becomes weaker 
and the exchange interactions J n decreases. 

Here we show that the calculated dependence of J„(f2) 
can explain the anti-Invar phenomenon in j-Fe^. If a 
magnetic subsystem is well described by the Heisenbcrg- 
like model {T]) its contribution to pressure according to 
the Hellman-Fcynman theorem is 



p — 



dHea 

on 



dJ, 



I -J 



an 



Zl m {e ° ei) 



(6) 



where Z\ is the number of nearest neighbors. We assume 
that the nearest-neighbor interaction is the strongest one 
which is supported by our fist-principle calculations. For 
a purpose of qualitative discussions, we will treat ex- 
change interactions perturbativcly assuming Z\J\ « T 



(T is the temperature and ks 
(e ei) w Ji/3T and therefore 



P„ 



7g dJf 



1). Then one has 



(7) 



This means that the pressure induced by magnetic ex- 
change interactions is positive and decreases with the 
temperature increase. 

The thermal expansion coefficient 



1 fan 



8P 



(8) 



can be divided into magnetic-structure independent part 
(ao) and one related with magnetic exchange interactions 
(a m ): a = ao + a m . The magnetic exchange part is equal 
to 



1 (dP„ 



B\dT 



B m ao 



(9) 



Here B is isothermal bulk modulus B = Bq + B m , Bq 
is the magnetic-structure independent part of B, and 
B m = —VL(dP m /dVl)T- Usually, the second term in 
Eq.© is neglected. Since dP m /dT < one can assume 
that the expression ([7]) should lead to the Invar behavior, 
and hence to the negative contribution to the thermal ex- 
pansion coefficient. A strong volume dependence of the 
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exchange parameter J\ can leads to the opposite conclu- 
sion. Substitutes Eq. ([7]) into Eq. ^ one finds 



1 Zi &7j 



1 - a Q T 



d\n{dJl/dQ) 



(10) 



Using calculated volume dependence Ji(fi) (Fig. [3]) and 
values ^ obtained from the experimental one can find 
that the second term in square brackets in the right-hand 
side of Eq. (Tit))) is approximately 1.1 at the temperature 
of 1200 K. Therefore, a total magnetic exchange contribu- 
tion to the thermal expansion coefficient ()10[) has positive 
sign. This corresponds to the anti-Invar behavior, in a 
qualitative agreement with the experimental data^l. 

The negative magnetic exchange contribution to the 
thermal expansion coefficient a m in the Invar materi- 
als usually is associated with a thermal dependence of 
the spontaneous magnetostriction^, while the positive 
contribution (anti-Invar behavior) is often considered to 
be related to thermal volume changes due to magnetic 
fluctuation o 24 i 55 . The present investigation allows us to 
explain the anti- Invar effect of high-temperature 7 phase 
of iron within a simple Heisenberg-like model in terms of 
magnetic softening of the bulk modulus, without any as- 
sumptions about two magnetic states of iron atoms with 
high and low volumes^. 

Due to the thermal expansion effective exchange pa- 
rameters increases with the temperature increase, 



J, 



ef 



Jo + AT, 



(11) 



formula into the mean-field expression for the magnetic 
susceptibility2i2£i2&, 



(12) 



3(T- 2J e/ /3) 



with a positive constant A > 0. If we substitute this 



one can see that corresponding temperature dependence 
leads to an increase of the effective magnetic moment, 
m 2 —¥ m 2 /(l — 2A/3) and the Curie temperature, Tq — > 
2b/(l - 2A/3). 

To conclude, we have carried out a systematic study of 
exchange parameters in a- and 7-Fe as functions of the 
volume and tetragonal deformation. The computational 
results demonstrate a strong coupling between lattice and 
magnetic degrees of freedom which should be taken into 
account in thermodynamic properties of Fe, especially its 
thermal expansion. Accurate analysis of the magnetic- 
structure independent contribution Eq allows us to con- 
clude that a response of fee and bec Fe to deformations 
is mainly controlled by the magnetic exchange. 
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